SU(3) versus deformed Hartree-Fock 
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Deformation is fundamental to understanding nuclear structure. We compare two ways to effi- 
ciently realize deformation for many-fermion wavefunctions, the leading SU(3) irrep and the angular- 
momentum projected Hartree-Fock state. In the absence of single-particle spin-orbit splitting the 
two are nearly identical. With realistic forces, however, the difference between the two is non-trivial, 
with the angular-momentum projected Hartree-Fock state better approximating an "exact" wave- 
function calculated in the fully interacting shell model. The difference is driven almost entirely by 
the single-particle spin-orbit splitting. 
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I. INTRODUCTION 

Deformation is one of the most salient features of 
atomic nuclei. The Bohr-Mottelson model, a quantum 
realization of the ellipsoidal surface of a liquid drop, was 
a major triumph in understanding nuclear structure jj]. 
The other major fundamental nuclear structure theory is, 
of course, the independent-particle or shell model, which 
is usually realized in a spherical basis. Elliott combined 
these two paradigms in the SU(3) model j|: the antisym- 
metric (fermionic) irreducible representations (irreps) of 
SU(3) in a spherical shell-model basis give rise to band 
structures associated with quadrupole deformed shapes. 

The nuclear Hamiltonian is not invariant under SU(3), 
however. Important parts of the nuclear Hamiltonian 
mix SU(3) irreps, like pairing and the single-particle spin- 
orbit splitting. To include symmetry breaking realis- 
tic interactions requires a more sophisticated version of 
the SU(3) model §, |, |]. The Nilsson model @, § of 
deformation describes single-particle levels by including 
quadrupole deformation and spin-orbit splitting, but one 
can think of Nilsson levels as a simplified version of de- 
formed mean-field, or Hartree-Fock, levels. 

It is commonly accepted that the simplest SU(3) the- 
ory, one that takes only the leading SU(3) irrep into ac- 
count, is a reasonable approximation to the exact ground 
state wavefunction of light sd-shell nuclei and not an un- 
reasonable first approximation as one moves to mid and 
upper s<i-shell nuclei and on into the lower /p-shell. This 
understanding provided motivation for a recent proposal 
for a "mixed-mode" approach to the selection of shell- 
model basis states, combining a few irreps of SU(3) to 
account for deformation with spherical configurations fa- 
vored by pairing-like modes and single-particle degrees 
of freedom, keeping the combined total number of basis 
states well-below that of the full space . 

While the mixed-mode approach has tremendous 



promise, SU(3) irreps are not the only means of realizing 
deformation. Projected Hartree-Fock states would be an 
obvious alternative. In fact, both SU(3) and projected 
Hartree-Fock have been made the basis of shell-model 
calculations ||, |j| However, a systematic compar- 
ison of SU(3) and Hartree-Fock as starting points for 
shell-model calculations has not beem made. This pa- 
per addresses that issue. Specifically, setting aside ques- 
tions related to the advantages or disadvantages of alge- 
braic versus numeric methods, we compare the leading 
SU(3) irrep and the angular-momentum-projected, de- 
formed Hartree-Fock state as approximations to "exact" 
shell-model ground state wavefunctions for sd- and lower 
/p-shell nuclei for "realistic" interactions. 

We find that projected Hartree-Fock is a better first 
approximation to the exact solution than the single lead- 
ing SU(3) irrep. Upon further investigation, however, we 
find that by removing the spin-orbit splitting of single- 
particle energies, the differences become almost negligi- 
ble, similar to the results found in [||. So although SU(3) 
has many advantages, being an algebraic theory, the re- 
sults of in this paper serve as a useful caveat. We also 
discuss some technical details of SU(3) versus projected 
Hartree-Fock, especially with regards to computational 
burden, in the concluding section. 



II. METHOD OF CALCULATION 

We start with a shell-model Hamiltonian which for the 
purposes of this paper we consider to be the true in- 
teraction. We compare the full numerical solution from 
shell- model diagonalization, \SM), against two approxi- 
mate wavefunctions: a state from the leading SU(3) irrep, 
\SU(3)), and the projected Hartree-Fock state \PHF). 
We consider only even-even nuclei. 



A. The model spaces 
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We work in complete Ohio model spaces, either the 
lsi/ 2 -0d3/2-0d 5 /2 (or sd) shell or the IP1/2-IP3/2-O/5/2- 
O/7/2 ipf) shell. This means there is an inert core, 16 



2 



for the sd-shell and 40 Ca for the pf shell. Valence parti- 
cles are restricted to the single-particle states of a single 
major harmonic oscillator shell. 



The interactions 



We use general one+two-body interactions which con- 
serve angular momentum and isospin. For the "true" in- 
teraction we use semi-empirical interactions such as the 
Brown- Wildenthal USD interaction for the sd shell fTlfl 
and the modified Kuo-Brown (KB3) interaction for the 
These interactions were derived by tak- 



pf shell (12 



ing a realistic G-matrix interaction computed from NN 
phase shifts, and then adjusting matrix elements to fit 
to a large number of levels over a wide range of A; for 
the USD interaction, this includes a phenomenological 
scaling of the two-body matrix elements by A~ - 3 (if one 
does not scale, the results below are qualitatively similar 
but quantitatively slightly weaker). What is important 
for our purposes is that these interactions describe very 
well the binding energies for a large number of nuclei and 
aside from the restriction to two-body, rotationally and 
isospin invariant interactions, were not restricted, that 
is, not assumed to have some schematic form such as 
quadrupole + pairing. Although less directly relevant 
to this paper, these interactions describe well excitation 
spectra and transition matrix elements. (Despite their 
general success, these general interactions are not perfect. 
For example, they have great difficulty in satisfactorily 
reproducing K-band splitting, something that is easily 
accomplished in an SU(3) basis by using a special scalar 
interaction that is of third+fourth order in the generators 
ofSU(3) @.) 

In the pf shell we also investigated the Brown- Richtcr 
interaction [0. The results were qualitatively the same, 
although the effect of spin-orbit splitting was smaller. 
For this paper we present only the KB3 results for the 
pf shell. 

Gueorguiev et al B projected exact shell-model wave- 
functions in the p/-shell onto SU(3) irreps and found that 
single-particle spin-orbit splitting made a significant con- 
tribution to mixing of SU(3) irreps in exact shell-model 
wavef unctions. Therefore for this study, analogous to ||, 
we removed the spin-orbit splitting of the single-particle 
energies for both the USD and KB3 interactions. 

Pairing also contributes to mixing of SU(3) irreps ||. 
Unlike single-particle spin-orbit splitting, however, pair- 
ing and Q ■ Q interactions share matrix elements in com- 
mon, so it is not trivial to disentangle pairing from Q ■ Q 
in realistic interactions such as USD and KB3. There- 
fore in this paper we focus only on the simpler case of 
spin-orbit splitting. 



C. "Exact" solutions 

The many-body Hamiltonians were constructed from 
the one+two-body interactions and diagonalized using a 
descendent of the Glasgow shell-model code [jl5| . 

In principle we allow all possible many-body configura- 
tions within that space; but because the Hamiltonian is 
rotationally invariant, we can without loss of generality 
restrict ourselves to those many-body states with M = 0, 
the so-called M-scheme basis. 



D. SU(3) solutions 

To construct the leading irreps of SU(3), we diago- 
nalized the SU(3) two-body Casimir C*2 — Q ■ Q + 3L 2 
(subtracting a small amount of L 2 in order to split states 
with the same (A, p) but different L). Because this is a 
two-body, real Hermitian operator, it can be treated as 
a Hamiltonian and therefore diagonalized by the same 
code that handles the realistic Hamiltonian. Unfortu- 
nately, because the shell-model code we used cannot han- 
dle three-body operators, we could not use the third- 
order SU(3) Casimir to split states with different (A, /i) 
but with degenerate values of Ci- Fortunately, for the 
leading irrep this was not a problem. 



E. Hartree-Fock solutions 

To construct the Hartree-Fock solutions, we took a 
Slater determinant wavefunction, and minimized the ex- 
pectation value of a given Hamiltonian using gradient de- 
scent. We then projected the Hartree-Fock Slater deter- 
minant onto the M = shell-model basis states used in 
computing the 'exact' and SU(3) solutions. This state did 
not have good angular momentum, but we projected out 
the J = part of the wavefunction. (We projected onto 
states of good angular momentum by a Lanczos moment 
method similar to that used in || to project onto SU(3) 
states and in [|l6] to generate Gamow- Teller strength dis- 
tributions. The unprojected state is used as the starting 
Lanczos vector, or "pivot," and then one peforms a few 
Lanzos iterations using J 2 as the Hamiltonian, splitting 
the pivot into components with good J. Although Ref. || 
used this technique to decompose shell-model wavefunc- 
tions among SU(3) irreps, for this paper we computed 
the wavefunction for the leading SU(3) irrep in the same 
shell model basis, as described in the previous subsection, 
and computed the overlap as a simple dot product, the 
wavefunctions being vectors in the M-scheme basis.) 

Our results for the ground state were insensitive to the 
orientation of the HF state. Results for excited states for 
triaxial nuclei were, however, sensitive to the orientation 
of the HF state. In general, results for excited states were 
qualitatively similar to the ground state, although quan- 
titatively less pronounced, so we do not present results 
for excited states here. 
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TABLE I: Energies of approximate wavefunctions, relative to 
the exact (shell-model) ground state, where HF = unprojected 
Hartree-Fock, PHF = projected Hartree-Fock, and SU(3) = 
leading irrep. All energies are in MeV. The label "nls" means 
single-particle spin-orbit splitting has been removed. 



Nucleus 


H 


AE HF 


AE PHF 


AE SU(3) 


20 Ne 


USD 


4.25 


0.95 


2.78 




USD-nls 


4.55 


1.02 


1.01 


24 Mg 


USD 


6.47 


1.75 


7.97 




USD-nls 


6.59 


1.45 


2.42 




USD 


7.16 


3.26 


10.07 




USD-nls 


6.12 


1.50 


2.36 


36 Ar 


USD 


3.95 


1.58 


5.29 




USD-nls 


3.70 


0.80 


0.84 


44 Ti 


KB3 


3.72 


1.18 


6.46 




KB3-nls 


3.80 


1.11 


2.54 



TABLE II: Overlaps of wavefunctions. SM= exact shell 
model ground state wavefunction, PHF = projected Hartree- 
Fock, and SU(3) = leading irrep. The symbol "nls" means 
single-particle spin-orbit splitting has been removed. 



Nucl 


H 


{SM\PHF) 


(SM\SU(S)) 


(SU(3)\PHF) 


20 Ne 


USD 


0.957 


0.881 


0.948 




USD-nls 


0.948 


0.952 


0.999 


24 Mg 


USD 


0.926 


0.667 


0.494 




USD-nls 


0.936 


0.892 


0.942 




USD 


0.847 


0.358 


0.171 




USD-nls 


0.925 


0.865 


0.944 


36 Ar 


USD 


0.923 


0.636 


0.586 




USD-nls 


0.948 


0.960 


0.973 


44 T j 


KB3 


0.933 


0.462 


0.473 




KB3-nls 


0.909 


0.727 


0.883 



confirmed by spectral decomposition as in g. This frag- 
III. RESULTS mentation is very sensitive to the single-particle spin- 

orbit splitting: a reduction of the spin-orbit splitting by 
A. Energies merely 15% will nearly triple (SU(3)\PHF). 



In Table | we compare the ground state energies, rela- 
tive to the exact (diagonalization of all valence configu- 
rations) ground state energy. That is, we compute 

AE = (*|-H"|*) - (SM\H\SM), 

where \SM) is the exact shell model wavefunction, and 
|^) = projected Hartree-Fock (PHF) or leading SU(3) ir- 
rep. We also show the energy of the unprojected Hartree- 
Fock (HF) for comparison. 

The expectation value of the exact Hamiltonian, with 
spin-orbit splitting, is considerably lower for the PHF 
state than for the SU(3) state, in one case by nearly 7 
MeV. When spin-orbit splitting is removed, the difference 
becomes much smaller, nearly vanishing for nuclei with 
axial symmetry. 



B. Wavefunctions 

In addition to computing the energies, we calculated 
the overlap of wavefunctions, shown in Table [n| The 
trends from the energies continue. For realistic inter- 
actions, the leading SU(3) has a smaller overlap with 
the exact shell model wavefunction than the projected 
Hartree-Fock state. When single-particle splitting is re- 
moved, the SU(3) state and the projected Hartree-Fock 
state come much closer to one another. The remaining 
difference may be due to pairing or other pieces of the 
realistic interaction. 

Note: the small overlap between the 32 S PHF state 
and the leading SU(3) irrep is due to the large fragmen- 
tation of the exact-shell model wavefunction over SU(3), 



C. Ground state geometry 

In addition to computing the overlap probabilities of 
wavefunction, it is useful to analyze the wavefunctions in 
terms of quadrupole deformation, i.e., the classic defor- 
mation parameters (3 and 7. 

We computed the deformation parameters /3, 7 for 
the unprojected Hartree-Fock state by diagonalizing the 
mass quadrupole tensor of the valence particles, as de- 
scribed by Ormand et al |P7| ]. For the interactions we 
chose the exact Hamiltonians (USD and KB3), the exact 
Hamiltonians with the single-particle spin-orbit splitting 
removed (USD-nls and KB3-nls), and the SU(3) second- 
order Casimir. (Note: although in principle one can com- 
pute deformation of SU(3) wavefunctions directly from 
(A, /i) j|, for consistency we also also applied Hartree- 
Fock to the SU(3) Hamiltonian.) 

Table HI shows that all three interactions yield similar 



results, but in particular the removal of spin-orbit split- 
ting makes the deformation geometry nearly identical to 
that of SU(3). This supports the overlap results from 
section HIB. We also see that inclusion of the spin-orbit 



force "softens" the deformation of the Hartree-Fock state, 
presumably through mixing of SU(3) irreps and similar 
to that seen in B. 



IV. DISCUSSION AND CONCLUSIONS 

In phenomenological but high-quality shell-model in- 
teractions, meaning that the one- and two-body matrix 
elements are fit to a large number of levels over a substan- 
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TABLE III: Deformation geometry of unprojected Hartree- 
Fock states for different Hamiltonians. (For SU(3) we used 
the second-order Casimir C2 as the Hamiltonian.) 



Nucleus 


H 





7 


20 Ne 


USD 


0.458 


0° 




USD-nls 


0.476 


0° 




SU(3) 


0.476 


0° 


24 Mg 


USD 


0.284 


13.7° 




USD-nls 


0.312 


18.6° 




SU(3) 


0.315 


19.1° 


32 g 


USD 


0.112 


24.6° 




USD-nls 


0.155 


40.4° 




SU(3) 


0.157 


40.9° 


36 Ar 


USD 


0.081 


60° 




USD-nls 


0.094 


60° 




SU(3) 


0.095 


60° 


44 Ti 


KB3 


0.440 


0° 




KB3-nls 


0.527 


0° 




SU(3) 


0.556 


0° 



tial range of A, a J-projected Hartree-Fock state does 
a nontrivially better job of approximating the "exact" 
ground state wavefunctions and binding energies derived 
from diagonalization in a full OHlo shell-model space, than 
does the leading SU(3) irrep. The difference is driven al- 
most entirely by the single-particle spin-orbit splitting, 
similar to that seen in previous work [Q pi. 

Does this mean that Hartree-Fock is necessarily supe- 
rior to SU(3)? The short answer is no; for a longer an- 
swer we compare HF and SU(3) along two axes: analytic 
insight and computational burden. 

SU(3) allows a tremendous degree of analytic insight 
into the structure of nuclear wavefunctions: for rotational 
nuclei both the ground state and excited states are often 
dominated by just a few SU(3) irreps, which can then 
illuminate the structure of band spectra including E2 
transition strengths. SU(3) is especially useful for sys- 
tems such as the rare earths [|l8|, where a full spherical 
shell-model calculation would be prohibitive. 

The computational burden can best be expressed 
by the number of basis Slater determinants needed to 
project out a good quantum number (such as J, M, A 
or fj.) , because modern shell-model codes use occupation- 
representation Slater determinants. In other words, one 
defines a single-particle basis and represents a many-body 
basis state as a binary word (l=occupied single particle 
state and 0=unoccupied). Much of the physics is con- 
tained by the choice of the single-particle basis. 

In a j-j coupled codes, such as Glasgow |H| , OXBASH 
HI, ANTOINE @, or REDSTICK /ELDORADO |T), 
the single-particle states couple I and s up to good j and 
to. An important advantage of j-j coupling is that spin- 
orbit splitting is trivial. The many-body basis states such 
single-particle states individually have good total M, and 



hence are called an M-scheme basis, but generally do not 
have good J. However one generally only needs a few 
hundred M-scheme Slater determinants to project out 
good J, as J 2 is block-diagonal in j-j basis. Deformation 
is harder, however: to project out good A or /i usually re- 
quires the entire many-body basis, which can be of order 
10 4 — 10 5 in the sd-shell and larger in the p/-shell. 

SU(3) shell model codes Q use L-S coupling built 
upon a cylindrical single-particle basis. In this basis spin- 
orbit splitting is hard (^requires many irreps) but de- 
formation is simple. For the leading irrep, to construct 
states of good A, \x one needs exactly one determinant, 
for the next-to-leading irrep the number is typically two 
or four depending upon whether the irrep occurs once 
or twice, and so on. Thus, because higher-weight irreps 
dominate the low-lying spectra, deformation is very sim- 
ple. In fact most of the computational burden lies again 
in projecting out states of good J, and is roughly similar 
to that for the spherical (j-j coupled) shell model. 

Finally we consider the Hartree-Fock state, which in 
principle incorporates both deformation and spin-orbit 
splitting. The HF state is not an exact eigenstate of 
SU(3) Casimirs, but that is not troubling as SU(3) is not 
an exact symmetry. Instead, as with other bases, the 
computational burden arises from projection of angular 
momentum. Formally one projects out states of good 
J by integrating over wavefunctions rotated by the Eu- 
ler angles weighted by Wigner D-functions. If, however, 
the HF state has axial symmetry (that is, good M) then 
one can project out a state of good J with just J max + 1 
Slater determinants, where J max is the largest total angu- 
lar momentum possible. For the sd-shell this is 14; for the 
p/-shell, this is 30. That is, one could get deformation, 
spin-orbit splitting, and rotational invariance with just a 
handful of Slater determinants. One could also imagine 
using particle-hole excitations on top of the Slater deter- 
minants and projecting onto good angular momentum, 
and indeed this approach has been used |9|, |Io[ ]. 

Unfortunately, for triaxial HF states one must in addi- 
tion project out states of good M; the number of states 
required goes from 10-30 to a few thousand. Further- 
more, HF states present additional difficulties when one 
considers multi-?iu; shells, where center-of-mass motion 
comes into play. Projection of a Hartree-Fock state onto 
a nonspurious state is difficult because the generators of 
center-of-mass motion are complicated and incomplete 
in a Hartree-Fock single-particle basis. The generaliza- 
tion of SU(3) to multi-fto; shells, the symplectic model 
p3| , automatically separates out spurious from nonspu- 
rious states. Therefore, while projected Hartree-Fock has 
some computational advantages, for multi-shell calcula- 
tions the symplectic extension of SU(3) may still be the 
method of choice. 
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